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Abstract 

i-^^l Results for the naive-time-reversal-odd quark distributions in a light-cone quark model are pre- 

I ' sented. The final-state interaction effects are generated via single-gluon exchange mechanism. The 

formalism of light-cone wave functions is used to derive general expressions in terms of overlap of 
. wave-function amplitudes describing the different orbital angular momentum components of the 

nucleon. In particular, the model predictions show a dominant contribution from S- and P-wave 
\ interference in the Sivers function and a significant contribution also from the interference of P and 

QQ ' D waves in the Boer-Mulders function. The favourable comparison with existing phenomenologi- 

0^ ■ cal parametrizations motivates further applications to describe azimuthal asymmetries in hadronic 

ly-^ . reactions. 
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I. INTRODUCTION 



Transverse momentum dependent parton distributions (TMDs), as an important exten- 
sion to the usual Feynman parton distributions, have attracted much attention in hadronic 
physics from both experiment and theory sides. Various hadronic processes have been used 
and proposed to study these distributions [l-21|. Together with the generalized parton dis- 



tributions (GPDs) (for reviews, see |22H26|). TMDs shall lead us to a comprehensive picture 
of parton distributions inside the nucleon, in particular, in a three-dimension fashion. 

Phenomenologically, in order to extract these distribution functions from experiments, 
we have to ensure that the QCD factorization applies in the associated processes. These 
issues have been extensively discussed in the last few years, and the relevant factorization 
theorem has been built up for a number of semi-inclusive processes, such as semi-inclusive 
hadron production in deep inelastic scattering and low transverse momentum Drell-Yan 



lepton pair production in hadronic collisions |27H29|. In the last few years, there has also 
been a remarkable experimental progress on experimental measurements (see Ref. [2l[ and 
references therein). More importantly, the proposed future experiments shall provide more 
constraints on these distribution functions. 

Meanwhile, reasonable model calculations of these transverse momentum dependent par- 



ton distributions have been proposed |3CH64| . These calculations promoted our understand- 
ing of the nucleon structure, and have been playing very important role as a first step to 
describe the experimental observations of the associated phenomena. In particular, these 
models provide us an intuitive way to connect the physical observables and the key input 
for the nucleon structure model, such as the quark spin and orbital angular momentum 
contributions to the proton spin. 

Transverse momentum dependent quark distributions are defined through the following 
quark-density matrix 

M{x,k^) = J ^^^^e-^'<{PS\m4^om\PS) , (1) 

where x and fc^ are the longitudinal momentum fraction and transverse momentum carried 
by the quark, respectively. Nucleon's momentum P is dominated by the plus component 
P"*" = (P° -|- P^) I \/2, and S represents the polarization vector. In the above equation, the 
gauge link L is very important to retain the gauge invariance and leading to nonzero naive- 
time-reversal-odd (T-odd) quark distributions. Among the leading order eight quark TMDs, 
six of them are called the naive-time-reversal even (T-even), whereas the rest two belong 
to the T-odd distributions. One is the so-called quark Sivers function, which describes the 
quark transverse momentum distribution correlated to the transverse polarization vector of 
the nucleon. The other is the so-called Boer-Mulders function, and usually interpreted as 
the transverse momentum correlated with the quark transverse polarization. Both quark 
distributions contribute to the azimuthal asymmetries in hadronic reaction processes. 



In Ref. , we have calculated the T-even quark distributions in a light-cone quark model 



extending previous works on the parton distribution functio ns ( PDFs) [35[, the GPDs 65 



681], nucleon form factors [69| and distribution amplitudes |70|. Such a model, based on 
the light-cone wave-function (LCWF) overlap representation, is able to capture the relevant 
information on the three-quark contribution to different observables. These calculations are 
well suited to illustrate the relevance of the different orbital angular momentum components 
of the nucleon wave function, and provide an intuitive picture for the physical meaning of 
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the quark TMDs. Moreover, they can be regarded as initial input for phenomenological 
studies for the semi-inclusive processes where quark TMDs play a very important role [7l|. 

In this paper, we extend these works to the T-odd quark distributions. The unique fea- 
ture for the latter distributions is the final/initial state interaction effects. Without these 
effects, the T-odd parton distributions would vanish. In the model calculation, these inter- 
actions are calculated by taking into account the one-gluon exchange mechanism between 
the struck quark and the nucleon spectators described by (real) LCWFs. This approach 



is complementary to a recent work [72j where the rescattering effects are incorporated in 
augmented LCWFs containing an imaginary phase which depends on the choice of advanced 
or retarded boundary condition for the gauge potential in the light-cone gauge. Recently, 
there has also been interesting study to g o beyond the one-gluon exchange approximation, 
by resumming all order contributions (6l|, 

The rest of the paper is organized as follows. In Sec. II, we briefly introduce the light- 
cone quark model, explaining its physical content and giving results for the light-cone wave- 
function amplitudes describing the different orbital angular momentum components of the 
nucleon state. In Sec. Ill, we derive the quark Sivers function. We present a general 
formalism in terms of overlap of light-cone wave-function amplitudes, and then apply it to 
a specific light-cone quark model wave function. The corresponding formalism for the Boer- 
Mulders function is described in Sec. IV. The model results for the T-odd distributions are 
presented in Sec. V and compared to different phenomenological parametrizations. Finally, 
we conclude with a section summarizing our findings. 



II. LIGHT-CONE AMPLITUDES IN A CONSTITUENT QUARK MODEL 

The wave-function amplitudes in light-cone quantization for the three-quark Fock state 
of the nucleon has been studied extensively in the literature According to the total 



quark orbital angular momentum projection, these wave-function amplitudes are classified 
into Iz = 0, Iz = 1, Iz = 2, /z = —1 components for total spin +1/2 of nucleon, i.e., 

I-P t)uud = \P t)uud + 1-^ t)w + 1-^ t)w + \P t)uud ■ (2) 

For completeness, we list the parametrization for these wave-function amplitudes following 



i'or complete] 
Refs. [ziS: 



1^ t)i=; = / d[l]dm3] ( 2, 3) + ^e^^h^h,iJ^l{l, 2, 3) 
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x^^It"(l) (^':(2)^t'(3) - ^V'(2)&l;(3)) |0) , (3) 



l^t)it = / d[l]dm3] [kt^i^^!,{l,2,3) + kt^^pl2AW) 



(&I;(l)^]r(2)^K(3) - bl,'{l)b]:i2)blli3)) |0) , (4) 



iPtY^:,-' = I d[l]d[2]d[3] fc-^V'il(l,2,3) 

x^6t«(i) {h]t{2)hll{?,) - 65:(2)4,(3)) |0) , (6) 

where a, /3 = 1,2 are transverse indexes and kf^^ = kf ± fcf. In Eqs. ([3])-(j6]) the integration 
measures are defined as 

^rii^roi^FQi dxidx2dx3 f \ d^ki^d^k2±d'^k3^ ( \ 

with Xi the fraction of the longitudinal nucleon momentum carried by the quarks, and 
ki±_ their transverse momenta. Furthermore, bl\ and b'^ ^ are creation and annihilation 
operators of a quark with fiavour q, helicity A and color i, respectively. In the following, 
we will describe the above light-cone wave-function amplitudes in a light-cone constituent 
quark model (CQM) following Ref. ji^. Working in the so-called "uds" basis 77, 78] the 
proton state is given in terms of a completely symmetrized wave function of the form 

l^' t) = \P t)w + \P t)udu + \P t)duu . (8) 

In this symmetrization, the state \P ']')udu is obtained from \P '\)uud by interchanging the 
second and third spin and space coordinates as well as the indicated quark type, with a 
similar interchange of the first and third coordinates for \P ^)duu- 

Following the derivation outlined in Ref. 65|, we find that the uud component of the 
light-cone state of the proton can be written as 

\P,^)uud = Y. I rf[l]^[2]rf[3]^;iik{^.fc.x;A.})^6U,(l)6i^^^ (9) 



x^&l;(l) (^5;(2)^I?(3) - (2)4^(3)) |0) , (5) 



In Eq. on]), assuming SU(6) spin-fiavor symmetry, we can factorize the LCWF '^uudii^i^ ^«-L! -^i}) 
in a momentum-dependent wave function and a spin-dependent part, i.e., 

Kud\{^^^ k^±■, A.}) = fc,^})i=|,^(Ai, A2, A3). (lO) 

In the above equation the momentum-dependent function is given by 



ip{{xi,k,^}) = 2{2n) 



1 U1U2UJ3 



1/2 



^Pi{x^,h^}), (11) 



Mo X1X2X3 

where ip{{xi, ki±}) is symmetric under exchange of the momenta of any quark pairs and 
is spherically symmetric, Ui is the free-quark energy, and Mq = Ui is the mass of the 
non-interacting three-quark system. The spin- dependent part in Eq. ( fTO!) is given by 

$a(Ai,A2,A3) = Yl (V2,/ii;l/2,/i2|l,/ii + /i2)(l,/ii + yU2;l/2,/i3|l/2,A) 

xDl{l{R,f{x,,k,^))Dl{ll{R,j{x2,k2^)^^^^^^ (12) 



In Eq. (fT2|) . D^^^[Rcf{x, k±)) is the matrix element of the Melosh rotation R^f 79 



{X\Rcf{x,kj_)\fi) 
m + xMq — icr 



(z X ki 



k- 



(13) 



The Melosh rotation corresponds to the unitary transformation which converts the instant- 
form spin eigenstates (the Pauli spinors) to light-front helicity eigenstates. In particular, 
the light-cone spin wave function of Eq. f|T2|) is obtained from the transformation of the 
canonical-spin wave function with zero orbital angular momentum component. The effects 
of the Melosh transformation are immediately evident in the presence of the spin-flip term 
z<T ■ (z X k±) in Eq. ( 1T3|) . Such a term generates non-zero orbital angular momentum, even if 
the original (instant-form) wave function only contained S-wave components. Therefore, as 
a consequence of total angular momentum conservation, the LCWF has components with 
total quark helicity different from the nucleon helicity. Making explicit the dependence on 
the quark helicities, the light-cone spin wave function of Eq. f lTS]) takes the following values: 



'i't(t,t,i)=n 

i 

^t(t4,t)=n 

i 

*ta,t,t) = n 

i 

(t, i, i) = n 

i 

*t (i, t, i) = n 

i 

*t(i,i,t)=n 

i 

^t(t,t,t)=n 

i 

(i,i,i)=n 
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1 


1 


^JN{Xi,ki 
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±) V6 
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±) V6 
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±) V6 
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^jN{xi, ki 


±) V6 



:(2aia2a3 + aik2 + (22^1 ^3")? 
— (ii(i20>^ ~h ^sk^ /c^ — '2,(iik^ k^ ^, 
(-010203 + a^kfk^ - 2a2kfk^), 
(0102^3" — ^r^^^a" ~ 2aia3/cJ), 
{—kfk2k^ + aia2k^ — 2a2a^ki), 



{a2a^ki + 0103/^^ + 2kfk2k' 
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-aiazk2 — a2Cizki +20102^:3), 
-a2kik^ — oi/cj/c;^ + 2a^kik2] 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



where Oj = (m -(- XjMo), and N{xi, ki±) = [{m + XiM^Y + 

Taking into account the quark-helicity dependence in Eqs. f|T4l) - fl2Tl) . the nucleon state 
can be mapped out into the different angular momentum components. After straightforward 
algebra, one finds the following representation for the nucleon wave-function amplitudes in 
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the light-cone CQM 

V^«(l,2,3) = ^({x,,fc,^}) 

X TT = — (-010203 + 03^1^ ■ k2± + 2aifci^ ■ k2± + 2aifc2x), 

(22) 

^(^)(1, 2, 3) = ^({x„ fc.x}) n , , 4(«3 + 2ai), (23) 

^/.(=')(1, 2, 3) = -^({x„ fc.x}) n , A («i«2 + fcL), (24) 

2, 3) = -^({x„ fc,x}) n , , 4^(QiQ2 + 2a3ai - fc^^ - 2fcix • fc2±), (25) 

i VA^(xi,fcix) v3 

^(^)(1, 2, 3) = ^^({x,, fc.x}) n , A («i«3), (26) 

^(^)(1, 2, 3) = V^({x„ fc.x}) n , , 4^2. (27) 

Notice that the results in Eqs. fl22p -f l27p follow from the spin and orbital angular momentum 
structure generated from the Melosh rotations, and are independent on the functional form 
of the momentum- dependent wave function. 

III. SIVERS FUNCTION 

The quark Sivers function can be calculated from the following definition 

fMx,kl) = -z{k^ + zky)^ I ^^^^e-^(«-^"-«-^-)(Ptl^(r,^J47^>Co^(o)l^4) • 

(28) 

As we discussed in the Introduction, the gauge link is crucial to obtain a non-zero Sivers 
function. In the covariant gauge, the gauge link can be reduced to the light-cone gauge 
linkEl. According to the light-cone wave function model, in the following calculations we 
choose the light-cone gauge A'^ = 0, where the gauge link reduces to a transverse gauge link 

at r 



AU+=o = ^exp \^-ig d'C± ■ Ax(r = 00, C±) j ■ (29) 

In the Sivers function of Eq. (128 p . we will expand the above gauge link to take into account 
the contribution from the one-gluon exchange diagram. Furthermore, in the light-cone gauge 
the gluon propagator takes the following form 

rf."(,) = -,- + !^fLip!. (30) 

[n-q] 



^ An ofF-light-cone gauge link has to be used to regulate the light-cone singularities for higher-order calcu- 
lations. In this paper, we will not encounter this singularity. Therefore, we will simply adopt the gauge 
link along the light-cone direction in covariant gauge and the transverse gauge link at spatial infinity in 
light-cone gauge. 
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FIG. 1. The leading contribution from the one-gluon exchange mechanism to the T-odd distribution 
functions. 



where n is the light-like vector ra^ = and n ■ q = . The gluon propagator has a light- 
cone singularity, as can be seen from the above equation. We will adopt the principal-value 
prescription to regulate this singularity. We have also checked that the final results do not 
depend on the prescription^ Under this prescription, there is no phase contribution from 
the above propagator. However, the transverse gauge link expansion, when combining with 
the n~q^/n ■ q factor of the above equation, leads to the following expression 

piq+oo 

= i7i6{q^) . (31) 



q+ 

This contribution provides the phase needed to generate a non-zero Sivers function. The 
dominance in the gluon propagator of the n~q^/n ■ q, with the _L index coming from the 
contraction with the transverse gauge link, also simplifies the interactions between the quark 
fields, since the quark scattering conserves the helicity. 

Finally, we obtain the following expression for the quark Sivers function 



f^-(x k') - _„2fclH^!*^ 1 1 /■ d*J£fc31 , , 



X 



E t l4^,(A:i)&,^,(A:)6l{3(A:3)6f,3(fc4)|Pi)}, (32) 

^-L Ai,A3 / i-j k,l 

where the quark momenta are defined a.s ki = k — q, = — q, T'* is the SUc{Q) Gell- 
Mann matrix and g is the gluon coupling with the quark field. Equation fl32l) corresponds 
to the diagrams in Fig. [1] with A = Ai and A4 = A3, for the helicity of the interacting and 
spectator quarks, respectively, and A = —A' for the helicity of the nucleon in the initial 
and final states. A few comments are in order to explain the above derivations. First, we 
have made an approximation for the interaction vertex between the gauge field from the 
gauge link and the quark fields in the proton wave function, by the covariant interaction 
form. In principle, we shall use the light-cone time-order perturbation theory to describe 
this interaction. However, we expect the modification being beyond the approximation we 
made in modelling the light-cone wave function itself. Nevertheless, it will be interesting to 
check how large these effects would be. Second, we used the perturbation theory to calculate 
the final-state interaction effects. For numerical estimate, we choose a reasonable value for 



^ For example, if we choose the so-called advanced boundary condition for the gauge potential, the transverse 
gauge link becomes unit, whereas the above gluon propagator generates phases which allow to recover the 
previous results with the principal-value prescription. 
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the strong coupling constant (see Sec. |V]). Meanwhile, we notice it may be not appropriate 
to use a perturbative coupling for this non-perturbative calculations. We regarded this as 
an important theoretical uncertainty, which exists in all model calculations of the Sivers 
function. 

As we discussed, in Eq. (132 p the quark helicity is conserved. On the other side, the 
hadron helicity flips from the initial to the final state. As a consequence, non-zero results 
for the Sivers function can be obtained only with a transfer of one unit of orbital angular 
momentum between the initial and the final nucleon states. 

Inserting in Eq. fl32|) the light-cone wave-function amplitude decomposition of the nucleon 
state introduced in Sec. [Tll one finds the following results in terms of the amplitudes 



/it ('^' ^- 



--gHd 
3^ 



P + zP fd^q^ 1 



(2vr)2 ql 



dx' / dH'^ / d[l]d[2\d[?>\^xiX2X^ 



(33) 



The function T^'^ for u quark is given by 



where 



^(1,2)^(3,4)^^^ 2, 3) - A(3'4)0(i.2)(i^ 2, 3) - A^^V^'Hl, 2, 3) + A(6)0(5)(l, 2, 3), (34) 



</.(i'2)(l, 2, 3) = ^«(1, 2, 3) - i{k\kl - k\kl)i,^''\\ 2, 3), 
0(3'^)(1, 2, 3) = fcrV^^'Hl, 2, 3) + k:,i,^'\\ 2, 3), 
0(^)(1, 2, 3) = kt^^''\\, 2, 3) - kt^^'^\\, 3, 2), 

<f\\, 2, 3) = fcr^3"V^^'Hl, 2, 3) - A;rA;2"^(^)(l, 3, 2). (35) 
The functions A in Eq. f lM|) are defined through 



^(1,2) 



b^ik - ki) 

r3 



5^{t' - k2)(j)^''^^*{2, 1", 3) + 5%t' - k3)(f)'^''''>*{2, 1", 3) 



(1,2)* 

+6\k - k2) y6\t' - ki) (^20(^-2)* (2", i, 3) + 0(^'''*(3, i, 2" 
+6^{t' - k^) (20(i'2)*(2", 1, 3) + </.(^'2)*(3, 1, 2")) 
+5\k - fcg) [d^'it' - k^) (^0(^'2)*(2, i, 3") + 0(i'2)*(3", i, 2) 
+6\t' - k^) f0(^'2)*(2, 1, 3") + 0(i'2)*(3", 1, 2) 



= 53^^ _ ^3(^/ _ ^^)0(3,4)*(2", 3) + ^3^^/ _ ^3)0(3,4)* (2", 1, 3) 

+5\k - ki) \5\t' - k2) f 2</.(=^'^)*(2, 1", 3) + 0(=''^)*(2, 3, l") 



+8\t' - fca) f2</.(=^'^)*(2, 1", 3) + </.(3'4)*(2, 3, l"))] 



+53(fc - k^) y5\t' - ki) (^</.(3'4)*(2, i, 3") + 0(3.4)*(2, 3", i; 

+5\t' - k2) f0(3'^)*(2, 1, 3") + 0(=^'^)*(2, 3", 1)^ 
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A(^) = 5''ik - ki) yS^'it' - k2) 2, 3) + 0(^)*(2, 1", 3) 

+5\t' - ks) 2, 3) + 0(^)*(2, 1", 3) 

+5'\k - k2) \5%t' - ki) f0(')*(i, 2", 3) + 



(2",i,3)) 



+5\t' - ks) 2", 3) + 0(5)* (2", 1, 3) 



= ^3^^ _ ^^3^^/ _ (^0(6)*(i"^ 2, 3) + 0(^)*(2, 1", 3) 
+6''{t' - k^) f0(6)*(i"^ 2, 3) + 0(6)* (2, 1", 3) 

■(i,2",3)- 



+5\k - k2) 5%t' - ki, 



(2",i,3)) 



+5\t' - ks 



■{1,2", 3) 



(2", 1,3) 



(36) 



where the quark coordinates are t" = {x,k± — qr_|_), and i = {x',t'j_ + S^{k — ki) = 
S{x — Xi)6'^{k± — ki±) and we used the notation 6^(t' — k^) = 6{x' — Xi)6'^{t'j_ — ki±). In the 
above equations, the complex conjugate only acts on the wave function ip^'^\ 

In Eq. (155]) . the contributions from the functions A^^'^^ and A^^'^) describe the interference 
between S and P waves, while the terms with A^^^ and A^^^ correspond to the contribution 
from P — D wave interference. 

Similarly for the d-quark, one has 

-p±d ^ 5(1,2)0(3,4)^^^ 2, 3) - 2, 3) - 5(5)0^*^^1, 2, 3) + 5(6)0(5)(1, 2, 3), (37) 

where the terms with i?(^'^) and 5(^,4) (^ggCTibg interference between S and P waves, while 
the terms with i?(5) a^fi ^(6) correspond to the contribution from P — D wave interference. 
The explicit expression for these functions is 

= 6^{k - ks) h^t' - A;2)0(^'2)*(2, 1, 3") + 6''{t' - k^)(f)'^^'^>{2, 1, 3") 
+S%k - fci) [S^{t' - k2) (0(''')*(2, 1", 3) + 0(^'2)*(3, 1", 2)) 
+6^{t' - ks) f0(''')*(2, 1", 3) + 0(^'')*(3, 1", 2) 



= s3(^k - ks) \s%t' - A;i)0(3'')*(2, i, 3") + S'\t' - A:2)0^'''^*(2, 1, 3") 
+6%k - k2) [S^if - fci) (0(3''^)*(2", i, 3) + 0(3'4)*(2", 3, i)' 
+6%t' - k^) f0(3'4)*(2", 1, 3) + 0(3.4)*(2", 3, 1)' 



5(5) = S%k - ks) ^6{t' - k2) (^0(5)*(1, 2, 3") + 0(5)*(2, 1, 3") 
+6{t' - ki) (0(5)* (i, 2, 3") + 0(5)* (2, i, 3") 



5(6) = 6^{k - ks) ^6^{t' - h) (0(')*(i, 2, 3") + 0(6)*(2, i, 3"^ 
+(5'(t' - A;2) f0(')*(l, 2, 3") + 0(6)*(2, 1, 3")^ 



(38) 



In the above equations, the complex conjugate only acts on the wave function ip^^\ 

Using the CQM expressions for the three-quark light cone amplitudes given in Sec. [TTl we 
obtain the following results for the Sivers function 



3^ fc? 



x6{x - X3)6\k^ - ks±)6{x' - xi)6\t'^ - k,^) r{{x[}, {k',^}) ^{{x,}, {h^}) 



dx' / dH'^ / d[l]d[2]d[3yxiX2X3 



X 



(39) 



where the quark momenta in the final state are (xg = x, k'^j_ = k^^ — q^] 



5 v-^i 



X ^ k 



t'j_ + {x2 = X2, k2j_ = k2±). In Eq. f l39l) . the functions and X^^ are given by 

3 

{k,}) = l[N~\k',)N~\h){zB,, + Bsy){A,A2 + B, ■ B2), 
1=1 
3 

X'\{k',},{k,}) = l[N~\k[)N~\h) 



(40) 



1=1 



l{- {A,A2 + B, ■ B2){iB,,, + Bsy) 

+ 2Bi ■ B3itB2, + B2y) + 2B2 ■ B-i{iBv, + B^y) 



+ 2i 



A3Ai{iB2^ + B2y) + ^3^2(^51^. + Biy) } 



(41) 



where 



Bi 



B. 



l,Z 



(m + x:M^)(m + XiMo) + k'^k^ + k!^k 
-(m + x[M'Q)k\ + (m + XiMo)k'y 
(m + x-Mo)fcf - (m + XiMo)fe 



Jx 



"i '"J '"I '"t 



(42) 



IV. BOER-MULDERS FUNCTION 

The calculation of Sec. IIIII can be repeated for the Boer- Mulders function, defined from 
the following quark correlation function 

hiix,kl)=^^ki^^ I ^^^e-^^^-''-^--'^-^lY.(P^\^^^^^^ ■ 

(43) 

Also in this case we expand the gauge link up to the next-to leading order, and following the 
same method we used in the calculation of the Sivers function, we find for the Boer-Mulders 
function 

h, ^(x, I ) = —q 7j — -— — = / — ^ = / d o 

' ^ ' kl 2 {2^Y^^(2k+)(2kt)J Vi^Hmt) J 



X 



^EEEE^F^'A.(i^A|6t^,(fcO&L(^)^IfA3(^3)&{.3(^4)|PA)}, (44) 



^-L A,A3 / i.j K 
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where the quark momenta are defined as ki = k — q, = — q. The above equation 
corresponds to the diagram of Fig. [1] with A = — Ai and A4 = A3 for the hehcity of the 
interacting and spectator quarks, respectively, and A = A' for the hehcity of the nucleon 
in the initial and final states, i.e. the helicity is conserved at the quark-gluon vertex, while 
the helicity of the struck quark fiips from the initial to the final state. Since the nucleon 
state has the same helicity in the initial and final state, the quark helicity fiip must be 
compensated by a transfer of one unit of orbital angular momentum. 

Inserting in Eq. dH]) the light-cone wave-function amplitude decomposition of the nucleon 
state introduced in Sec. [Til one finds the following results in terms of the amplitudes 



{2ny ql 



d[l\d[2\d[^yxiX2X^'H- 



(45) 



where the function 'H-^'^ for the up quark is 

= _c(i'2)0(3.4)(i^ 2, 3) + 2, 3) - 2, 3) 

+C^(6) 0(3,4)^^^ 2, 3) + 2, 3) - 2, 3), 



(46) 



with 



2, 3) = V^(i)(l, 2, 3) + i{klkl - k\kl)^^^\\, 2, 3), 
0(=^'^)(1, 2, 3) = k\^^'\\, 2, 3) + kti^^'\\, 2, 3), 
0('^)(1, 2, 3) = ktkt^^''\\, 2, 3) - A;+A;2+^(6)(l, 3, 2). 



(47) 



In Eq. f l46|) . the terms containing C*-^'^-* and C^^'^ describe the contribution from S and 
P wave interference, while C-^-* and C*-^-* are associated with the P — D wave interference 
terms. The explicit expression for these functions is 

C7(i'2) = ^3^^ _ J^3(^/ _ (^0(i.2)*(i^ 3^ 2") + 20(i'2)*(2", 3, i)) 
+6^{t' - ks) (0(''')*(1, 3, 2") + 20(i'2)*(2", 3, 1))" 
+6%k - ks) \6''{t' - ki)(f)^^'^>{3", 2, i) + 6%t' - A;2)0(^'')*(3", 2, 1) 



Ci3A) = s3^k - ki) ^6^{t' - k2) (0(=''')*(3, 2, 1") + 20(3,4)*^3^ ^//^ 2)) 
+ 6\t' - ks) (0('''^*(3, 2, 1") + 20(=^'')*(3, 1", 2))" 
+6^{k - ks) \6^{t' - ki)4>^^'^>{l, 3", 2) + 6%t' - A;2)0^''^^*(i, 3", 2) 

C(3.4) = s^^k - ki) [5\t' - A;2)0(''")*(l", 2, 3) + 5\t' - A;3)0(''")*(l", 2, 3)] 
+5\k - k2) [S%t' - A;i)0(=^'^)*(2", i, 3) + S'it' - k^)^^'''> {2" , 1, 3)] 

(7(6) =S^{k- ki) ^6''{t' - k2) (0(*^)*(1", 2, 3) + 0(*^)*(2, 1", 3)) 
+ 5'(t' - ^3) f0^'^*(l", 2, 3) + 0(^)*(2, 1", 3)' 
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= S%k - h) \s%t' - /t2)0^'''^*(2, 1", 3) + S%t' - h)4>^''^>i2, 1", 3) 
+6\k - k2) W'it' - ki)4>^^'^>{l 2", 3) + ^^(t' - k3)4>^^'^^*{l, 2", 3) 



C;(5) = s%k - k2) [S%t' - ki) 2", 3) + 0(5)*(2", i, 3) 

+6''{t' - ks) f0(^)*(l, 2", 3) + 0(5)*(2", 1, 3) 



(48) 



In the above equations, the complex conjugate only acts on the wave function 
Analogously, the function T-L^ for the down quark is 

^ /}(i-2)0(3,4)^^^2,3) -£)(3,4)^(i,2)^^^2,3) + ^(3,4)^(6)^^^2,3) 
_^(6) 0(3,4) 2, 3) + 2, 3) - D^'^4>^''^\1, 2, 3), 



(49) 



where the S — P wave interference contribution comes from the terms proportional to Z^^^'^^ 
and D^^'^\ while the remaining two terms give the contribution from the interference of P 
and D waves. The function D in Eq. ( l49l) are defined as 

D(i.2) = ^3^^ _ ^^3(^/ _ A;i)0(i'2)*(i^ 2, 3") + 53(t' - A;2)0(^''^*(l, 2, 3")] , 



^(3,4) ^ ^3^^ _ r^3(^/ _ A;,)0(3'4)*(i^ 2, 3") + 6'' {f - /ts)^^''^)* (1, 2, 3") 

= ^3^^ _ [<53(t'-A;i) (</.(=^'^>(3",2,i) + 0(3,4)*(3'/J^2)) 
+ S'{t'-k2) (0(3'^)*(3",2,l) + 0(3,4)*(3'/^i^2))] , 



= s^^k - ki) (t' _ (0(')*(3, 2, 1") + 0(^)*(2, 3, 1"; 
+ (5=^(t' - ks) f0('^*(3, 2, 1") + 0(6)* (2, 3, 1"; 



= 6%k - k2) \5\t' - k^) (0(^'')*(i, 2", 3) + 0(i'2)*(3, 2", i)) 
+ (^'(i' - fca) f0^'''^*(l, 2", 3) + 0(^'2)*(3, 2", 1)' 



= 53^^ _ J^3(^/ _ (^0(5)*(i, 2", 3) + 0(^)*(3, 2", i) 

+5\t' - k^) (0(^)*(1, 2", 3) + 0(5)*(3, 2", 1))] . (50) 

In the above equations, the complex conjugate only acts on the wave function ip^^\ 

In the model for the three-quark light cone amplitudes introduced in Sec. Ull we find the 
following explicit results 



P - iky f d^q^ 1 



dx' / dH'^ / d[l]d[2]d[3yxiX2X3 



x6{x - X3)S\k^ - ks^)6{x' - xi)5\t '^-k r{{x-}, {K^}) {k,^}) 
X35.3., |5,,i/2X°0({fc:},{fc,}) + ^[5.^1/2 + 25.,-i/2]Xi^({fc:},{fc.})|, (51) 
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where the quark momenta in the final state are (xg = x, k'^j_ = fcsj, — qj_] 



(4 



X2,k2j_ = k2±). In Eq. ( 151]) . the functions X^^ and are given by, 



't=i 

3 



X'\{k'},{k}))=llN~\{k'm^\{k,}) 



1=1 



X - I {3A,A2 - B, ■ B2) A3 + 2{A,B2,. 



A2B,JB: 



3,x 



+2{A,B2,y + A2B,^,)B3,y + 2(Ai52,, + A2B,^,)B3,, 

where the functions Ai and Bi are defined in Eq. f l42|) . and 

is = (m + X3Mo)(fc^" + 2^^) - (m + x'3M^)(A;^ + ik^), 



-i{m + x'^M'^){m + X3M0) + + + ifcf), 

(m + x'3M^)(m + X3M0) + (4,, + z/c^,)(fc^ + zA;|), 
2(m + x;5M^)(A:^ + ikl) + i(m + x^Mo){k'^ + i^^). 



(52) 



(53) 



(54) 



V. RESULTS AND DISCUSSION 



The formahsm described in the previous sections is apphed in the following to a specific 
CQM, adopting a power-law form for the momentum-dependent part of the light-cone wave 
function, i.e. 

N' 

^({x.fc.x})= ^^2^^,^^ , (55) 

with A^' a normalization factor. In Eq. f l55|) . the scale /3, the parameter 7 for the power-law 
behaviour, and the quark mass m are taken from Ref. |80| . i.e.. /3 = 0.607 GeV, 7 = 3.4 
and m = 0.267 GeV. According to the analysis of Ref. [Sli these values lead to a very 
good description of many baryonic properties. The same parametrization of the momentum 
dependent part of the LCWF in Eq. flS^ has been successfully a ppl ied also in recent works 
for the calculation of the electroweak properties of the nucleon |69l |, 
and T-even TMDs |4i,[7l|. 

In order to fix the coupling constant appearing in Eqs. (l33i) and (H5 



GPDs 26, 35, 65-67 



we need to determine 

the hadronic scale of the model. This is achieved in a model independent way following the 
prescription of Ref. [s^], by matching the value of the momentum fraction carried by the 
valence quarks, as computed in the model, with that obtained evolving backward the value 
experimentally determined at large Q^. The strong coupling constant 05 (<5^) entering the 
evolution code at NLO is computed by solving the NLO transcendental equation numerically. 



In 



A2 



An 



(3! 



In 



A' 

(31 



0, 



(56) 



as obtained from the renormalization group analysis [82|, |83 
expression 



It differs from the more familiar 



1 



Att PoHQV^ 



NLO J 



^ (3'o 



f3i lnln(gVA 



NLOy 



ln(QVA 



NLO J 



(57) 
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FIG. 2. Results for the first transverse-momentum moment of the Sivers function, for up (left) 
and down (right) quark, as function of x. The dashed curves show the results at the hadronic 
scale of the model /Iq = 0.094 GeV'^, and the solid curves correspond to the results after NLO 
evolution to = 2.5 GeV^, using the evolution pattern of the unpolarized parton distribution. 
The lighter and darker shaded areas are the uncertainty bands due to the statistical error of the 
parametrizations of Ref. [50(] and Ref. [30|, ISJ], respectively. Both parametrizations refer to an 



average scale of = 2.5 GeV^ 



valid only in the limit Q"^ ^ 



^^■NLO' 



where Anlo is the so-called QCD scale parameter. 



The hadronic scale, /Xq, consistent with the presence of valence degrees of freedom only 
is /iQ = 0.094 GeV^, with Anlo = 0.248 GeV. This corresponds to a value of the strong 
coupling constant in Eq. fISB]) as{fJ'l)/{4:'7i) = g'^/{4:nY = 0.14, and is consistent with the 
analysis of Refs. |46l - where a similar procedure was adopted. 

The first transverse-momentum moments of the Sivers and Boer-Mulders functions are 
shown in Figs. [2] and El using the definition 



k- 



2M2 



(58) 



with j = f^j? and j = hf, respectively. In the figures the dashed curves correspond to the 
results at the hadronic scale of the model /Xq, while the solid curves are obtained by applying 
a NLO evolution to = 2.5 GeV^, assuming for the first transverse-momentum moment 
of the Sivers function the same anomalous dimension of the unpolarized parton distribution 
and for the first transverse-momentum moment of the Boer-Mulders the evolution pattern 
of the chiral-odd transversity distribution. Although these are not the exact evolution 
patterns, this is the standard procedure adopted so far in model calculations |45l-l48| and 
parametrizations 30 . Isol. of the T-odd TMDs, since the exact evolution equations are 



still under study [l7, 84-88] and evolution codes for these distributions are not yet available. 
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For the Sivers function in Fig. [2] we also show the results from recent parametrizations, 
valid at an average scale of = 2.5 GeV^, obtained from a fit to available experimental 
data on transverse single spin asymmetries for pion and kaon in semi-inclusive deep inelastic 
scattering. In particular, the darker shaded area represents the uncertainty due to the 
statistical errors in t he p arametrization of Ref. (Hof, while the lighter shaded area corresponds 
to the same for Ref. The model predictions for the contribution of u and d quarks are 



of the same order of magnitude and opposite sign, and after evolution are well compatible 
with the phenomenological parametrizations. The effects of the evolution are crucial to 
reproduce the position of the peak at x ~ 0.2 for both the u and d quark distributions, and 
to rescale the magnitude of the distributions within the range of the parametrizations. 

A non trivial constraint in model calculations of the Sivers function is given by the 
Burkardt sum rule 8^ 



E 



q=u, d,s,g,- 



dxft/^\x,kl) = 0, (59) 



which corresponds to require that the net (summed over all partons) transverse momentum 
due to final-state interaction is zero Restricting the sum in Eq. flS^ over the up- and 



down-quark contributions, our model calculation of the Sivers function reproduces exactly 
the sum rule. 

In Fig.[3]we compare the model results for the absolute value of the Boer-Mulders function 
with phenomenological parametrizations obtained from recent fits to available experimental 



data. In particular, the dashed-dotted curve corresponds to the analysis of Refs. [91|, |92| at 
the average scale of = 2.4 GeV^ of the HERMES Q and COMPASS B lo^ measure- 
ments of the cos 20 asymmetry in SIDIS, while the short-dashed curve shows the results of 
Refs. 9^, 97 1 valid at ~ 1 GeV^, obtained from a fit to pd j98| and pp 9^ Drell-Yan data 



measured by the E866/NuSea Collaboration, with the shaded area describing the variation 
ranges allowed by positivity bounds. We note that the available data do not allow yet a 
full fit of hj; with its x and k]_ dependence and these phenomenological parametrizations 
are only first attempts to extract information on this distribution. Upcoming experimental 
SIDIS data also from JLab and plans for Drell-Yan experiments at GSI will play a cru- 
cial role to better constrain these analysis. Our model predictions after the "approximate" 
evolution to = 2.4 GeV^ are compatible with the phenomenological analysis of SIDIS 
data, reproducing both the peak position and the behaviour in x, while are at variance with 
the analysis of the Drell-Yan data. In particular we confirm the find i ngs of Ref. [9l| and 
the expectations from various theoretical analysis j4o|, 45, 48, 59, lOd - 102 |. predicting the 



same sign for both the up and down contributions, with the u component of hj; larger in 
magnitude than the corresponding component of and the d components of and 
with approximately the same magnitude and opposite sign. 

In Fig. m we show the decomposition of the Sivers and Boer- Mulders functions in the 
contributions from the different partial-wave amplitudes of the nucleon LCWF. The dashed 
curves correspond to the results from the interference of S and P waves, the dotted curves 
show the contribution from P — D wave interference, and the solid curves are the total 
results, sum of all the partial wave contributions. The S — P wave interference terms give 
the dominant contribution to the Sivers function of both u and d quarks, while the P — D 
wave interference terms contribute at most by 20% of the total results. On the other side, 
the relative weight of the P — D wave interference terms increases in the case of the Boer- 
Mulders function. It corresponds to 30% of the total results for the up-quark distribution 
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FIG. 3. Results for the first transverse-momentum moment of the Boer-Mulders function, for 
up (left) and down (right) quark, as function of x. The dashed curves show the results at the 
hadronic scale of the model /Xg = 0.094 GeV^ and the solid curves correspond to the results after 
NLO evolution to = 2.4 GeV^, using the evolution pattern of the transversity distribution. 
The dashed-dotted curves are the results of the phenomenological parametrization of Refs. |9ll. |92| 
at the average scale of = 2.4 GeV^, and the short-dashed curves correspond the results of 
Refs. 1961,1971] valid at ~ i GeV^, with the shaded area describing the variation ranges allowed 



by positivity bounds. 



and becomes the dominant contribution in the case of down quark, reac hing up to 60% of the 
total resuh. We also note that, contrary to the case of T-even TMDs [4^, the assumption 
of SU(6) symmetry in the model does not imply any proportionality between the T-odd 
distributions of up and down quark. As outlined in Ref. j48|, this is due to the fact that 
in the case of the T-odd functions one is using a two-body operator associated with FSI, 
while for the T-even TMDs the proportionality results from the calculation with a one-body 
operator. In comparison with other model calculations, our light-cone model predictions are 
similar in shape but significantly different in magnitude from those in the non-relativistic 



CQM of Refs. |46H48| . The main differences in this calculation can be traced back to the use 
of covariant quantization and non-relativistic wave functions. Furthermore, the quark-gluon 
interaction vertex is treated non relativistically. Analogous discrepancies are evident in the 



comparison of our predictions with the results of the bag-model j47|, |48|, l58| , although in this 



case the calculation is fully relativistic. Since the bag model uses covariant quantization, 
the spin structure is worked out in terms of canonical spin instead of light-cone helicity, 
and at the quark-gluon vertex one has contributions from both spin-flip and spin-conserving 
terms. However, these terms reduce to helicity conserving terms (A3 = A4 in the diagram of 
Fig. Hj) when written in terms of light-cone helicity, in agreement with our model calculation. 
For a more detailed discussion on the relation between the structure of TMDs in terms of 
canonical spin and light-cone helicity we refer to Ref. 
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0.2 0.4 0.6 0.8 



FIG. 4. Angular momentum decomposition of the first k± moment of the Sivers function for the 
up (left panel ) and down (right panel) quark. The dashed curves show the contribution from 
the interference of S and P waves, and the dotted curves correspond to the contribution from the 
interference of P and D waves. The solid curves are the total results, sum of all the partial- wave 
contributions. 



Finally, with respect to the diquark models of Refs. |40|, |45| we have different magnitude and 
shape for both the Sivers and Boer Mulders functions. The different magnitude might be 
due to the choice of different values for the quark-gluon coupling constant as well as to the 
absence of D-wave components in these calculations. Note however that our results are at 
variance with the calculations in the diquark models also for the relative magnitude between 
up- and down-quark contributions. 

The dependence on x and kl_ of the Sivers and Boer-Mulders functions is shown in Figs. E] 
and El respectively, for the separate up (left) and down (right) quark contributions. The 
behaviour in is very similar for the two distributions, and does not depend on the quark 
flavour. The fc^-dependence shown in Figs. [5] andj6] is definitely not of Gaussian form. 
However, following the exercise performed in Ref. [7l| for the T-even distributions, it is 
interesting to compare the model predictions for the mean square transverse momenta with 
the results of the Gaussian model. We define the mean transverse momenta {n = 1) and 
the mean square transverse momenta (n = 2) for the TMD j{x, k]_) as follows 

_ Jdxjd'k^ klj{x,kl) 



J dx J d'^k± j{x,k^) 
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FIG. 5. The Sivers function as function of x and k\ for up (left panel) and down quark (right 
panel) . 




where k± = The corresponding results for the T-odd distributions are shown in Table I. 
In the Gaussian model the following relation holds 

{kl) ^^^^ - {k^)^ , (61) 
vr 

which implies that the ratio shown in the last column of Table I should be equal to one. The 
model results deviates from unit by 10%. We also note that the mean transverse momenta 
in Table I are quite small, much smaller than expected from phenomenological studies. This 
is due to the low scale of the model, and Sudakov effects are expected to make the fc^ 
distributions larger when evolving to larger and experimentally relevant scale. 

In Fig. [71 we show the spin density in the transverse-momentum space of unpolarized up 
(left panel) and down (right panel) quark in a transversely polarized nucleon, defined as 



(62) 



with S± the nucleon transverse-polarization vector, and fi{x, kjj the monopole distribution 
corresponding to spin densities for unpolarized quarks in an unpolarized target. When S± 
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TMD j 


{k± 


) in GeV 




in GeV^ 


TTikl) 


frr 
hi 


up 
0.22 
0.23 


down 
0.24 
0.24 


up 
0.071 
0.077 


down 
0.084 
0.080 


up 
0.90 
0.90 




down 
0.90 
0.91 



TABLE I. The mean transverse momenta and the mean square transverse momenta of T-odd 
TMDs, as defined in Eq. ()60p . from the hght-cone CQM. If the transverse momenta in the TMDs 
were Gaussian, then the result for the ratio in the third column would be unity, see text. 



points in the x direction, the dipole contribution related to the Sivers function introduces a 
large distortion on the monopole term, perpendicular to both the spin and the momentum 
of the proton and with opposite sign for up and down quarks. The corresponding average 
transverse-momentum shift is defined as 



y\i 

J IT 



f d^k^kyp"^^ {k 



(63) 



and results 



M 

IT 2 



-70.31 MeV, 



JIT 



M 



dxf[^- 



(x) 



140.62 MeV. 

(64) 



The fact that the absolute value of the average transverse momentum induced by the 
Sivers function is twice as large for d quark than for u quark is just a consequence of 
the Burkardt sum rule in Eq. (l59ll . This intrinsic k± shift is the an alog ous of the dipole 



10J|. Although it is 



deformation related to the GPD E in impact-parameter space [63|, 
not possible to establish a model independent relation between the Sivers function and the 
GPD E [63[, we note that the LCWF overlap representation of E, for vanishing longitudinal 
momentum transfer, is given in terms of the same combinations of light-cone amplitudes 
parametrizing the Sivers function in the one gluon-exchange approximation, but evaluated 
for different values of quark variables (65, The values for the average shifts in impact- 
parameter space within the present light-cone quark model were found (6^)" = k"/(2M) = 
0.20 fm and {b^)'^ = k'^/{M) = — 0.33fm j68[, where k'^ is the quark contribution to the 
proton anomalous magnetic moment. 



Analogously, the spin density of transversely polarized quarks and unpolarized nucleon 
is related to the Boer-Mulders effect by 



dx- 
2 



ffix^kD + sV^k^—hl^ix^k 



M 



(65) 



where s± is the quark transverse-polarization vector. In Fig. [8] we show the spin density 
for quark polarization in the x direction. Since the Boer-Mulders function is negative for 
both up and down quarks, the sideway shift is always in the positive y direction. The 
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FIG. 7. Spin density in the transverse-momentum plane for unpolarized quarks in a transversely 
polarized nucleon. The left panel is for up quark, and the right panel for down quark. 
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FIG. 8. Spin density in the transverse-momentum plane for transversely polarized quarks in an 
unpolarized nucleon. The left panel is for up quark, and the right panel for down quark. 



corresponding average dipole distortion is 

M 

T 



159.40 MeV, {ky)t = M 



-215.73 MeV. 

(66) 



Although the Boer-Mulders function is smaller in magnitude for down quark than for up 
quark, one observes that the average sideways distortion for down quark is stronger. This is 
because the monopole distribution related to ff is twice as large for up quarks as for down 
quarks, therefore adding the dipole contributions results in a more pronounced distortion 
for down quarks than for up quarks. The corresponding dipole distribution in impact- 



parameter space is described by the chiral odd GPDs Et + 2Ht- As found in Ref. (67 
these GPDs for zero longitudinal momentum transfer are given by the same combination 
of LCWFs which enter the calculation of in the one-gluon exchange approximation, but 
at different kinematics. The corresponding average distortion in impact-parameter space is 
proportional to the tensor anomalous magnetic moment nip, and in the present light-cone 
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quark model is g iven by k^/(2M) = 0.42 fm and k,^/{M) = 0.55 fm for up and down quark, 
respectively |68j. 

VI. CONCLUSIONS 

In this paper we have investigated the naive-time-reversal-odd quark distributions, the 
quark Sivers and Boer-Mulders functions, in a light-cone quark model. The final-state in- 
teraction effects are calculated by approximating the gauge link operator with a one-gluon 
exchange interaction. In this framework, we have derived the general formalism for the T-odd 
quark distributions in terms of overlap of light-cone wave function amplitudes describing the 
different orbital angular momentum components of the nucleon state. This model indepen- 
dent expressions are particularly suitable to emphasize the correlations of quark transverse 
momentum and transverse polarizations of the nucleon and of the quark. For numerical 
estimates, the nucleon light-cone wave-function has been constructed by assuming a light- 
cone constituent quark model with SU(6) spin-flavor symmetry and a momentum-dependent 
part which is spherically symmetric. Under this assumption the orbital angular momen- 
tum content of the wave function is fully generated by the Melosh rotations which boost 
the rest-frame spin into the light-cone. As a result, we found explicit expressions for the 
light-cone amplitudes which match the analytic structure expected from model-independent 
arguments 74-7^1. The model dependence enters the choice of the momentum-dependent 



part of the light-cone wave function. In this work, we adopted a phenomenological descrip- 
tion, by assuming a specific functional form with parameters fitted to hadronic structure 
constants. The same wave function was used to predict many other hadronic properties, 
providing a good description of available experimental data, and being able to capture the 
main features of the quark contribution to hadronic structure functions, like parton distribu- 
tions jsll, generalized parton distributions j65-68||, nucleon form factors 69 1 and distribution 



amplitudes j70| , and T-even transverse momentum dependent quark distributions ji^, 7ll . 
The corresponding results for the Sivers and Boer-Mulders function have been presented in 
this paper by showing the decomposition into the contributions from different orbital angu- 
lar momentum components. Both functions require a transfer of one unit of orbital angular 
momentum between the initial and final states. In particular, the Sivers function for both 
up and down quark is dominated by the interference of S- and P-wave components, while 
the P — D wave interference terms contribute at most by 20% of the total results. On the 
other side, the relative weight of the P — D wave interference terms increases in the case of 
the Boer-Mulders function, in particular for the down-quark component. Furthermore, the 
model results for the Sivers function satisfy exactly the so-called Burkardt sum rule, which 
is a non-trivial constraint for model calculations and parametrizations. 

In order to compare with phenomenological parametrizations obtained from a fit to avail- 
able experimental data for semi-inclusive deep inelastic scattering and Drell-Yan processes, 
we evolved the model results to the experimental scale. Since the exact evolution equations 
for the T-odd quark distributions are still under study, we used those evolution equations 
which seem most promising to be able to simulate the correct evolution. We evolved the 
first transverse-momentum moment of the Sivers function by means of the evolution pattern 
of the unpolarized parton distribution, while for the first transverse-momentum moment of 
the Boer-Mulders we used the evolution pattern of the transversity. After evolution, the 
model results are consistent with the available parametrizations, especially for the Sivers 
function. There is agreement between the signs of the various flavor components, and also 
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for the magnitude and the position of the maxima in x. These findings encourage further 
phenomenological apphcations of the model to describe azimuthal asymmetries in hadronic 
reactions. 

We also found that the x and k\ dependence is similar for the Sivers and Boer-Mulders 
functions, and approximately independent on the quark flavor. In particular, the k\ is 
not of Gaussian form. However, it is worthwhile to evaluate the degree of approximation 
introduced by the Gaussian Ansatz within the model in the calculation of observables. This 
task is left for future applications of the model. 

Finally, we discussed the spin densities in the transverse-momentum space related to the 
Sivers and Boer-Mulders effects, showing that they are consistent with the model results 
for the corresponding spin densities in the impact-parameter space described by generalized 
parton distributions. 
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